Rounding of a first-order quantum phase transition to a strong-coupling critical point 
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We investigate the effects of quenched disorder on first-order quantum phase transitions on the 
example of the Af-color quantum Ashkin-Teller model. By means of a strong-disorder renormaliza- 
tion group, we demonstrate that quenched disorder rounds the first-order quantum phase transition 
to a continuous one for both weak and strong coupling between the colors. In the strong coupling 
case, we find a distinct type of infinite-randomness critical point characterized by additional in- 
ternal degrees of freedom. We investigate its critical properties in detail, and we discuss broader 
implications for the fate of first-order quantum phase transitions in disordered systems. 
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The effects of disorder on quantum phase transitions 
have gained increasing attention recently, in particular 
since experiments have discovered several of the exotic 
phenomena predicted by theory (see, e.g, Refs. [1, 2]). 
Most of the existing work has focused on continuous tran- 
sitions while first-order quantum phase transitions have 
received less attention. 

In contrast, the influence of randomness on pure 
systems undergoing a classical first-order transition 
has been comprehensively studied. Using a beauti- 
ful heuristic argument, Imry and Wortis [3] reasoned 
that quenched disorder should round classical first-order 
phase transitions in sufficiently low dimension and thus 
produce new continuous phase transitions. This analy- 
sis was extended by Hui and Berker [4]. Aizenman and 
Wehr [5] rigorously proved that first-order phase transi- 
tions cannot exist in disordered systems in dimensions 
d < 2. If the randomness breaks a continuous symmetry, 
the marginal dimension is d = 4. 

The question of whether or not disorder can round a 
first-order quantum phase transition (QPT) to a contin- 
uous one was asked by Scnthil and Majumdar [6], and, 
more recently, by Goswami et al. [7]. Using a strong- 
disorder renormalization group (SDRG) technique, they 
found that the transitions in the random quantum Potts 
and clock chains [6] are governed by the well-known 
infinite-randomness critical point (IRCP) of the random 
transverse-field Ising chain [8, 9]. The same holds for the 
A^-color quantum Ashkin-Teller (AT) model in the weak- 
coupling (weak interaction between the colors) regime [7] . 
This implies that disorder can indeed round first-order 
quantum phase transitions. 

In the strong-coupling regime of the AT model, on the 
other hand, the renormalization-group (RG) analysis of 
Ref. [7] breaks down. Goswami et al. speculated that 
this implies persistence of the first-order QPT in the 
presence of disorder, requiring important modifications 
of the Aizenman- Wehr theorem. However, shortly af- 
ter, Greenblatt et al. [10, 11] proved rigorously that the 
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FIG. 1. (Color online) Schematic of the renormalization- 
group flow diagram in the disorder — coupling strength pa- 
rameter space. For e < e c (left arrows), the critical flow ap- 
proachs the usual Ising infinite-randomness critical point of 
Ref. [9]. For e > e c (right arrows), we find a distinct infinite- 
randomness critical point with even stronger thermodynamic 
singularities. 



Aizenman- Wehr theorem also holds for quantum systems 
at zero temperature. 

In this Letter, we resolve the apparent contradiction 
between these results. We show that quenched disor- 
der rounds the first-order QPT of the AT model in the 
strong-coupling regime as well as in the weak-coupling 
regime. Moreover, we unveil a distinct type of infinite- 
randomness critical point governing the transition in 
the strong-coupling regime. It is characterized by ad- 
ditional internal degrees of freedom which appear be- 
cause a higher symmetry is dynamically generated at 
criticality. As a consequence, the critical point dis- 
plays even stronger thermodynamic singularities than the 
transverse-field Ising IRCP. To obtain these results, we 
have developed an implementation of the SDRG method 
that works for both weak and strong coupling. In partic- 
ular, this method can deal with the diverging intercolor 
interactions as well as the associated additional degen- 
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eracies. A schematic of the resulting RG flow is shown 
in Fig. 1. 

In the remainder of this Letter, we sketch the deriva- 
tion of our renormalization-group method, we determine 
its fixed points, and we contrast the resulting critical be- 
havior in the weak and strong-coupling regimes. 

The Hamiltonian of the one-dimensional TV-color quan- 
tum AT model [12-14] is given by 
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Here, i indexes the lattice sites, a and {3 index colors, 
and cr x and a z are the usual Pauli matrices. The interac- 
tions Ji and transverse fields hi arc independent random 
variables taken from distributions restricted to positive 
values, while th,i and eji (also restricted to be positive) 
parametrize the strength of coupling between the colors 
[15]. Various versions of the AT model have been used to 
describe layers of atoms absorbed on surfaces, organic 
magnets, current loops in high-T c superconductors as 
well as the elastic response of DNA molecules. Note the 
invariance of the Hamiltonian under the following dual- 
ity transformation: cr z al a z al+1 -> T x i7 o x a i -> r z iT z i+1 , 
Ji hi, and ejj <^ eh,%, where t x and t z are the dual 
Pauli operators. In the clean version of our system with 
N > 3, the QPT between the paramagnetic and ordered 
(Baxter) phases is of first-order type [12-14, 16]. 

In order to tackle the Hamiltonian (1), we now develop 
a SDRG method. In the weak-coupling regime (eh, £j "C 
e c , where e c is some TV-dependent threshold), our method 
agrees with that of Goswami et al. [7] . Here, we focus on 
the strong coupling regime e/j, ej S> e c where the method 
of Ref. [7] breaks down. 

The basic idea of the SDRG method consists in iden- 
tifying the largest local energy scale and perturbatively 
integrating out the corresponding high-energy degree of 
freedom. As we are in the strong-coupling regime, this 
largest local energy is either a four-spin interaction ( "AT 
interaction") fcj = tj t iJi or a two-color field-like term 
("AT field") Qi = th.ihi- We thus define our high-energy 
cutoff il = maxjfci, g{\. 

We now derive the decimation procedure. If the largest 
local energy is an AT field located, say, at site 2, the 
unperturbed Hamiltonian for the decimation of this site 
reads 
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The ground state (GS) of H is \(f>o) = ■ ■ ■ 

with energy Eq = —Ng 2 /eh, 2 — N(N — 1)52/2, where 
each arrow represents a different color. Flipping n colors 
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FIG. 2. (Color online) Spectrum of the unperturbed Hamil- 
tonian (2) as function of the number of colors flipped with 
respect to the ground state |— ¥,—¥, •■■ ,—*■)■ As long as 
T < fl/eh, the pseudo ground state \<f>' ) = \<— ,•<— ) 

can be neglected when computing observables (stage 1 of the 
RG). When T > Q/th, |<^o) and \<j>' ) become effectively degen- 
erate implying that both states need to be taken into account 
(stage 2 of the RG). 



leads to (^) degenerate excited states with energy E n = 
Eq + 2ng 2 /eh,2 + 2n(N — n)g 2 . In the strong-coupling 
regime, eh S> 1, the state \<f/ ) = \<— ■ ,<—) plays a 
special role. Its energy E' = Ng 2 /ch, 2 — N(N — 1)32/2 
differs from that of the true ground state only by the 
subleading Ising term E' — E = 2Ng 2 /ef l _ 2 (see Fig. 2). 
It can thus be considered a "pseudo ground state" which 
may be important for a correct description of the low- 
energy physics. The true and pseudo ground states each 
have their own sets of low-energy excitations which we 
call the ground-state and pseudo-ground-state sectors of 
low energy states. 

The couplings of site 2 to its neighbors, 
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is the perturbation part of the Hamiltonian. We now dec- 
imate site 2 in the 2nd order perturbation theory, keeping 
both the true ground state and the pseudo ground state. 
It is important to note that 2nd order perturbation the- 
ory does not mix states from the two sectors as long as 
N > A [17]. After decimating site 2, the effective inter- 
action Hamiltonian of the neighboring sites reads (in the 
largc-e,/ limit) 
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FIG. 3. (Color online) a) Decimating a site results in a renor- 
malized bond (characterized by k and ej) between its neigh- 
bors, and it introduces an extra binary "sector" degree of 
freedom represented as an Ising spin ( = ±1 in an external 
field uj [see Eq. (7)]. b) Decimating a bond results in a renor- 
malized site characterized by g, eh and another sector degree 
of freedom. 



Here, ( = ±1 is a new Ising degree of freedom which 
represents the energy splitting between the true and the 
pseudo ground states. In the large-ej regime, it is only 
very weakly coupled to the rest of the chain and can be 
considered free. In Fig. 3(a), we sketch this decimation 
procedure. 

The decimation of a bond can be treated in the same 
way. If an AT four-spin interaction, say fo, is the largest 
local energy, the unperturbed Hamiltonian reads 
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Its GS is obtained by any sequence of parallel nearest 
neighbors pairs (e.g. \<j> ) = | tt, ft, 44, 44= tt. 44, • • • , tt 
)) with energy E = -TVfc 2 /ej )2 - TV(TV - l)fc 2 /2. As 
above, in the strong-coupling limit ej$ 3> 1, Hq has a 
pseudo-GS consisting of a sequence of anti-parallel nearst 
neighbors pairs (e.g. \</>' ) = | tt, U, It, It, tt, tt, • ■ ■ , tt 
)) with enegy E' = Nk 2 /ej - TV (TV - l)fe/2. 

When integerating out the bond, the two-site cluster 
gets replaced by a single site which contains one addi- 
tional internal binary degree of freedom, namely, whether 
the cluster is in the GS sector or in the pseudo-GS sector. 
Its effective Hamiltonian reads 
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Here, ( distinguishes the two sectors as before. The du- 
ality of the Hamiltonian can be seen by comparing Eqs. 
(5) and (8) after exchanging the roles of k and g as well 
as eh and ej. 

Note that the magnetic moment ft of the new effec- 
tive site depends on the internal degree of freedom £ [see 
Fig. 3(b)] because neighboring spins are parallel in the 



GS sector but antiparallel in the pseudo-GS sector. We 
will come back to this point when discussing observables. 

The SDRG proceeds by iterating these decimations. In 
this process, the coupling strengths ej, eh flow to infinity 
if their initial values are greater than some e c (TV). This 
means that the Ising terms J; , hi become less and less 
important with decreasing energy scale f2. The large-e 
approximation thus becomes asymptotically exact. The 
remaining energies are the AT four-spin interactions ki 
and the AT fields <?i. Their recursions relations have the 
same multiplicative structure as the recursions of Fisher's 
solution [9] of the random transverse-field Ising model. 
The flow of the distributions P(ki), R{gi) and their fixed 
points are thus identical to those of Fisher's solution, 
see Fig. 1. We conclude that the distributions of k, g 
have an infinite-randomness critical fixed point featuring 
exponential instead of power-law scaling [8, 9, 18]. As 
the Ising coupling Ji, hi have vanished, this critical fixed 
point has the symmetry of the AT interaction and field 
terms which is higher than that of the full Hamiltonian. 

The zero-temperature phase diagram of our system is 
determined by the low-energy limit of the SDRG flow. 
There are three classes of fixed points parameterized by 
the distance from criticality r = (In g) — (In k) (where 
(• • • ) denotes the disorder average): The critical fixed 
point at r = 0, and two lines of fixed points for the or- 
dered (r < 0) and for the disordered (r > 0) Griffiths 
phases. This implies that there is a direct continuous 
phase transition between the ordered (Baxter) and disor- 
dered phases. We found no evidence for additional phases 
or phase transitions. In agreement with the Aizenman- 
Wehr theorem [10], we thus conclude that disorder turns 
the clean first-order QPT into a continuous QPT in both 
strong-coupling and weak-coupling regimes. 

We now turn to the behavior of observables at low 
temperatures. The basic idea is as follows [19]. We dec- 
imate the system until the cutoff energy scale Q reaches 
the temperature T. For low enough T, the distributions 
of all energy scales in the renormalized system become 
very broad, and thus, the remaining degrees of freedom 
can be considered as free. Applying this procedure, we 
have to distinguish two stages depending on the impor- 
tance of the pseudo ground state. (1) Both AT and Ising 
couplings are above the temperature. (In this stage, we 
decimate sites and bonds whose internal sector degrees of 
freedom are frozen in the true ground state, C = 1-) (2) 
The temperature is below the AT couplings but above the 
Ising couplings. (Here we still decimate sites and bonds, 
but their internal degrees of freedom are free, i.e., they 
can be in either of the two sectors, £ = ±1.) 

Let us illuminate this RG scheme on the example of 
the entropy [20]. When the RG flow stops at fi = T, 
all spins are completely free. A surviving cluster has 
2 N available states (2 per independent color) giving an 
entropy contribution of TV In 2, i.e., Schain = wtTVui2, 
where no is the density of surviving clusters at energy 
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scale [9]. Moreover, during stage 2 of the flow, resid- 
ual entropy was accumulated in the internal degrees 
of freedom, each of them contributing In 2 to the en- 
tropy. Noticing that each stage-2 RG decimation gen- 
erates one extra degree of freedom, and that stage 2 
starts when fi/ej^ = T, (e x is the typical value of e x ^ 
at energy scale fi), the extra contribution to the entropy 
is Septra = [wj(n ejT - n T ) + w h (n ehT - n T )] In 2, with 
wj = 1 — Wh being the fraction of coupling decimations 
in the entire stage 2 of the RG flow. To compute 5 e xtra 
we need to know how ej and e/j depend on Q. From the 
recursions (8) and (5), it is clear that lne/j (and lnej) 
scale like the number of sites (bonds) in a renormalized 
cluster (larger bond). 

At criticality, wj = w h = 1/2, n n ~ [^(Oj/O)]- 1 /^, 
with "0 = 1/2 being the tunneling exponent, and lne^ = 
lnej - pn(n//n)]*, with <j> = i(l + VE) [9, 21]. Thus, 
summing the two contributions we find that 



S = d 



T 



(9) 

where C\ and C2 are nonuniversal constants, and Vti is 
the bare energy cutoff. As > 1, the low-T entropy 
becomes dominated by the extra degrees of freedom S — > 

Sextra-KOz/T)]" 1 /^. 

In the ordered Griffiths phase (r < 0), wj — > 1 and 
lnej = A^W-^lntfij/fi), with A being a nonuni- 
versal constant of order unit, v = 2 the correlation 
length exponent, and z = l/(2|r|) the dynamical expo- 
nent. As no ~ \r\ v ip,/^li) 1 / z , we find that SWtra ~ 
\r\ v (T /VLi) 1 /( z+Az ' 1 '} ln2, which dominates over the chain 
contribution ~ T X / Z N In 2. As expected from duality, the 
same result holds for the disordered phase (r > 0). 

To discuss the magnetic susceptibility, we need to find 
the effective magnetic moment fj, e s of a cluster surviv- 
ing at the RG energy scale Q, = T . If all internal de- 
grees of freedom were in their ground state, jj, e s would 
be given by the number of sites in the cluster. How- 
ever, analogously to the entropy, ^t e ff is modified because 
of the stage 2 of the RG flow. In this stage, the inter- 
nal degrees of freedom are free, meaning not all spins 
in a surviving cluster are parallel, reducing the effective 
moment. A detailed analysis based on the central limit 
theorem [22] gives /i e ff ~ [hi(f2//T')] <i ' / ' 2 at criticality and 
u eff ~ r^-^HOj/T)] 1 / 2 in the disordered Griffiths 
phase, as well as fj, e g ~ r - ^/ 2 T -1 /( 2z ) in the ordered 
Griffiths phase. 

The magnetic susceptibility x(T) can now be com- 
puted. All eliminated clusters had AT fields greater 
than the temperature, and thus do not contribute to 
X since they are fully polarized in cc-direction, whereas 
the surviving clusters are effectively free and contribute 
with a Curie term: x ~ ^eS^/T. We find that 
X ~ [ln(0 / /T)] < ^- 1 /^/T in the critical region, while it be- 
comes x ~ r u+2u ^ 1 -^T 1 / z ~ 1 ln(fi//T) in the disordered 



Griffiths phase, and take a Curie form x ~ \r\ v ~^T~ l in 
the ordered Griffiths phase. 

In summary, we have solved the random quantum 
Ashkin- Teller model by means of a strong-disorder 
renormalization-group method that works not just for 
weak-coupling but also in the strong-coupling regime and 
yields asymptotically exact results. In the concluding 
paragraphs, we put our results into broader perspective. 

First, we have demonstrated that random disorder 
turns the first-order QPT between the paramagnetic and 
Baxter phases into a continuous one not just in the 
weak-coupling regime but also in the strong-coupling 
regime. This resolves the seeming contradiction between 
the quantum Aizenman-Wehr theorem [10, 11] and the 
conclusion that the first-order transition may persist for 
sufficiently large coupling strength [7] . 

The resulting continuous transition is controlled by two 
different IRCPs in the weak and strong coupling regimes. 
For weak coupling, the critical point is in the universality 
class of the random transverse- field Ising chain [9]. For 
strong coupling, we find a distinct type of IRCP which 
features a higher symmetry than the underlying Hamilto- 
nian. The associated internal degrees of freedom lead to 
even stronger thermodynamic singularities both at criti- 
cality and in the Griffiths phases. 

Our results apply to N > 4 colors where the true and 
pseudo ground-state sectors are not coupled. As a result, 
the Ising terms in the Hamiltonian are irrelevant pertur- 
bations (in the renormalization group sense) at our IRCP. 
The case N < 4 is special because the two sectors get 
coupled and thus requires a seperate investigation. In- 
terestingly, novel behavior has been recently verified for 
the classical transition in the two-dimensional AT model 
[23] ioT N = 3. 

We note that our method is also interesting from a gen- 
eral renormalization-group point of view. After a deci- 
mation, the resulting system cannot be represented solely 
in terms of a renormalized quantum AT Hamiltonian be- 
cause the internal degree of freedom needs to be taken 
into account. Normally, the appearance of new variables 
dooms an RG scheme [24]. Here, however, the new vari- 
ables, despite their influence to observables, are inert in 
the sense that they do not influence the RG flow of the 
other terms in the Hamiltonian, which makes the problem 
tractable. We expect that this insight may be applicable 
to renormalization-group schemes in other fields. 

The strong-disorder RG approach to the random quan- 
tum AT-modcl gives asymptotically exact results for both 
sufficiently weak and sufficiently strong coupling (c<l, 
e 3> 1), see Fig. 1. The behavior for moderate e is 
not exactly solved. In the simplest scenario, the weak- 
coupling and strong-coupling IRCPs are separated by a 
unique multicritical point at some intermediate coupling, 
however, more complicated scenarios cannot be excluded. 
The resolution of this question will likely come from nu- 
merical implementations of the SDRG and /or (quan- 
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turn) Monte-Carlo simulations. 
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